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Abstract. In cores of supernovae and crusts of neutron stars, nuclei can adopt interesting shapes, such as rods or slabs, etc., 
which are referred to as nuclear “pasta.” Recently, we have been studying the pasta phases focusing on their dynamical aspects 
with quantum molecular dynamic (QMD) approach. We review our findings on the following topics: dynamical formation 
of the pasta phases by cooling down the hot uniform nuclear matter; a phase diagram on the density versus temperature 
plane; structural transitions between the pasta phases induced by compression and their mechanism. Properties of the nuclear 
interaction used in our works are also discussed. 


I. INTRODUCTION 

In ordinary matter, atomic nuclei are roughly spherical. This may be understood in the liquid drop picture of the 
nucleus as being a result of the forces due to the surface tension of nuclear matter, which favors a spherical nucleus, 
being greater than those due to the electrical repulsion between protons, which tends to make the nucleus deform. 
When the density of matter approaches that of atomic nuclei, i.e., the normal nuclear density po, nuclei are closely 
packed and the effect of the electrostatic energy becomes comparable to that of the surface energy. Consequently, 
at subnuclear densities around p ~ Po/2, the energetically favorable configuration is expected to have remarkable 
structures: the nuclear matter region (i.e., the liquid phase) is divided into periodically arranged parts of rodlike or 
slablike shape, embedded in the gas phase and in a roughly uniform electron gas. Besides, there can be phases in 
which nuclei are turned inside out, with cylindrical or spherical bubbles of the gas phase in the liquid phase. These 
phases with nonspherical nuclei are often referred to as nuclear “pasta” phases because nuclear slabs and rods look like 
“lasagna” and “spaghetti.” Likewise, spherical nuclei and bubbles are called “meatballs” and “cheese,” respectively. 

In equilibrium dense matter in supernova cores and neutron stars, existence of the pasta phases has been predicted by 
Ravenhall et al. JH and Hashimoto el al. El . Since these seminal works, properties of the pasta phases in equilibrium 
states have been investigated with various nuclear models. They include studies on phase diagrams at zero temperature 
El Iz5l l33l l3hl I 4 H 1 and at finite temperatures El . These earlier works have confirmed that, for various nuclear 
models, the nuclear shape changes as: sphere —> cylinder —> slab —> cylindrical hole —> spherical hole —> uniform, with 
increasing density. 

In these earlier works, however, a liquid drop model or the Thomas-Fermi approximation is used with an assumption 
on the nuclear shape (except for Ref. I 4 H 1 L Thus the phase diagram at subnuclear densities and the existence of 
the pasta phases should be examined without assuming the nuclear shape. It is also noted that at temperatures of 
several MeV, which are relevant to the collapsing cores, effects of thermal fluctuations on the nucleon distribution 
are significant. However, these thermal fluctuations cannot be described properly by mean-field theories such as the 
Thomas-Fermi approximation used in the previous work El- 

In contrast to the equilibrium properties, dynamical or non-equilibrium aspects of the pasta phases had not been 
studied until recently except for some limited cases 114, 26]. Thus it had been unclear even whether or not the pasta 
phases can be formed and the transitions between them can be realized during the collapse of stars and the cooling of 
neutron stars, which have finite time scales. 

To solve the above problems, molecular dynamic approaches for nucleon many-body systems are suitable. They 
treat the motion of the nucleonic degrees of freedom and can describe thermal fluctuations and many-body correlations 


beyond the mean-field level. 

Using the framework of QMD Qj], which is one of the molecular dynamic methods, we have solved the following 
two major questions 0ES1I. 

• Question 1: Whether or not the pasta phases are formed by cooling down hot uniform nuclear matter in a finite 
time scale much smaller than that of the neutron star cooling? 

• Question 2: Whether or not transitions between the pasta phases can occur by the compression during the collapse 
of a star? 

The pasta phases have recently begun to attract the attention of many researchers (see, e.g.. Refs. nm and refer¬ 
ences therein). The mechanism of the collapse-driven supernova explosion has been a central mystery in astrophysics 
for almost half a century (e.g.. Ref. 01). Previous studies sugg est that the revival of the shock wave by neutrino heating 
is a crucial process. As has been pointed out in Refs. j3^|38] and elaborated in Refs. 1 12;(3S3, the existence of the 
pasta phases instead of uniform nuclear matter increases the neutrino opacity of matter in the inner core significantly 
due to the neutrino coherent scattering by nuclei 00 ; this affects the total energy transferred to the shocked matter. 
Thus the pasta phases could play an important role in the future study of supernova explosions. Our recent work |37| 
strongly suggests the possibility of dynamical formation of the pasta phases from a crystalline lattice of spherical 
nuclei; effects of the pasta phases on the supernova explosions should be seriously discussed in the near future. 


II. METHOD: QUANTUM MOLECULAR DYNAMICS 

Among various versions of the molecular dynamic models, Quantum Molecular Dynamics (QMD) (TJ is the most 
practical one for investigating the pasta phases. Rodlike and slablike nuclei are mesoscopic entities of nuclei them¬ 
selves and they contain a large number of nucleons. QMD, which is a less elaborate in the treatment of the exchange 
effect, allows us to study such large systems with several nonspherical nuclei. The typical length scale r c of half of the 
inter-structure is r c ~ 10 fm and the density region of interest is around half of the normal nuclear density p (l . The total 
nucleon number N necessary to reproduce n structures in the simulation box is N ~ po(2r c n) 3 (for slabs). It is thus 
desirable to use ~ 10 4 nucleons in order to reduce boundary effects. Such large systems are difficult to be handled by 
other molecular dynamic models such as FMD @] and AMD 12411 . whose calculation costs increase as ~ A 4 , but are 
tractable for QMD, whose calculation costs increase as ~ N 2 . 

It is also noted that the exchange effect is less important for the nuclear pasta structures, which are in the 
macroscopic scale for nucleons. This can be seen by comparing the typical values of the exchange energy and of 
the energy difference between pasta phases. Suppose there are two identical nucleons, i = 1 and 2, bound in dif¬ 
ferent nuclei. The exchange energy between these particles is calculated as an exchange integral: I = f U(r\ — 
Vi) <pi (ri)<p^(r 2 )<P 2 (ro)<p >2 ( r l) dridrj , where U is the potential energy. An asymptotic form of the wave function 
is given by <p,- ~ exp {—hr) with kt = QlmEi/h , (i = 1,2), where is the binding energy and m is the nucleon 
mass. The exchange integral reads 1 ~ exp [— {k\ + U)/'’] ~5x 10 6 MeV for the internuclear distance R ~ 10 fm and 
Ei ~ 8 MeV, which is extremely smaller than the typical energy difference per nucleon between different pasta phases 
of order 0.1 keV (for neutron star matter) - 10 keV (for supernova matter). Therefore, it is expected that QMD is not a 
bad approximation for investigating the pasta phases. 


1. Model Hamiltonian and its Properties 

In our studies on the pasta phases, we have used a nuclear force given by a QMD model Hamiltonian with the 
medium-equation-of-state parameter set in Ref. El . This model Hamiltonian consists of six parts: 

— K -\- Vpauli V Vskyrme Vsym V V \] |) -f Vcoulomb , (1) 

where K is the kinetic energy; Vp au ii is the momentum-dependent “Pauli potential,” which reproduces the effects of the 
Pauli principle phenomenologically; Vskyrme is the Skyrme potential which consists of an attractive two-body term and 
a repulsive three-body term; V sym is the symmetry potential; Vmd is the momentum-dependent potential introduced as 
two Fock terms of the Yukawa interaction; Vcoulomb is the Coulomb energy between protons. 

The parameters in the Pauli potential are determined to fit the kinetic energy of the free Fermi gas at zero 
temperature. The above model Hamiltonian reproduces the binding energy of symmetric nuclear matter, 16 MeV 



per nucleon, at the normal nuclear density po = 0.165 fm 3 and other saturation properties: the incompressibility is 
set to be 280 MeV and the symmetry energy is 34.6 MeV. This model Hamiltonian also well reproduce the properties 
of stable nuclei relevant to our interest: the binding energy except for light nuclei from 12 C to 20 Ne [lSj], and the 
rms radius of the ground state of heavy ones with A > 100 115]. It is also confirmed that another QMD Hamiltonian 
close to this model provides a good description of nuclear reactions including the low energy region (several MeV per 
nucleon) (23il . 

Let us then examine other properties of the nuclear interaction (at zero temperature), which have not been determined 
accurately yet but have important effects on inhomogeneous structure of matter at subnuclear densities. Such quantities 
are the nuclear surface tension £ sur f, the energy per nucleon e n of the pure neutron matter, and the proton chemical 
potential pp 0> in the pure neutron matter. The surface tension E sur f, which is the most important among these three 
quantities, controls the size of the nuclei and bubbles, and hence the sum of the Coulomb and surface energies. With 
increasing £ sur f and so this energy sum, the density p m at which matter becomes uniform is lowered. There is a 
tendency, especially in a case of neutron star matter, that the higher e n , p m is lowered. This is because larger e n tends 
to favor uniform nuclear matter without dripped neutron gas regions than mixed phases. In neutron star matter, there is 
also a tendency that the lower Pp°*, the smaller p m . This is because - Up 0 represents the degree to which the neutron 
gas outside the nuclei favors the presence of protons in itself. 



FIGURE 1. (Color) The nuclear surface energy per unit area (the surface tension) E sur r versus the proton fraction X[ n in the 
nuclear matter region. The red solid squares are the values of the present QMD Hamiltonian [I19i]; the solid curve is the result of the 
Skyrme-Hartree-Fock calculation with a modified version of 1’ parameter set done by Ravenhall, Bennett and Pethick (RBP) IH; 
the dotted curve is from Baym, Bethe and Pethick (BBP) 0. This figure is adapted from Ref. I. 3ljl . 

In Figs. [fland|3 we have plotted the results of these quantities for the present model Hamiltonian. We can say that, on 
the whole, they give reasonable values within uncertainties of each quantities. It is noted that E sm f of the present model 
shows moderate values between the RBP and BBP results. The behavior of Up’ 1 is also in a reasonable agreement with 
various Skyrme-Hartree-Fock calculations except for higher densities of p n > 0.1 fm~ 3 relevant only for neutron star 
matter just below p m . The quantity e n , however, shows some deviation from the major trend of Skyrme-Hartree-Fock 
results and of the values of the microscopic calculations. The behavior of e n of the present model is similar to that of 
the SkM interaction and of the interaction by Myers et al. [22]. 

In the following, for each case of different value of the proton fraction x of matter, we summarize the consequences 
of the features of Esurf, £« and of the present QMD interaction. 

1. For symmetric nuclear matter (x — -Vn — 0.5) 

According to E sur f at jq n = 0.5, the present model is consistent with the other results, and is an appropriate effective 
interaction for studying the pasta phases at x = 0.5. 
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FIGURE 2. (Color) The neutron density p n dependence of the energy per nucleon e n (left panel) and the proton chemical 

pote ntial Pp°^ (right panel) of pure neutron matter. The red solid squares show the result of the present QMD model Hamiltonian 
QjJ- The red dotted line denoted by SLy4 is the result from Ref. Q], and the colored broken lines as marked by the other Skyrme 
interactions (FPS21, T, FPS and SkM) are the results summarized by Pethick, Ravenhall and Lorenz m . The black dotted line is 
the result of Sjoberg 01 and the black dashed line is that of Myers, Swiatecki and Wang SI . The open squares shows the result 
of the GFMC calculation by Carlson et al. (g|; the open stars denote the values obtained by Akmal, Pandharipande and Ravenhall 
d; the triangles are those from Friedman and Pandharipande |l(J; the crosses from Siemens and Pandharipande 13 11 1. The large 

error bars and the scatter of Pp 0 ^ of QMD in the low density region p < 0.3po are due to the local rou ghne ss of the density in the 
neutron matter, which is caused by the fixed width of the wave packet. This figure is adapted from Ref. 1331 . 


2. For neutron star matter (x < 0.1) 

The melting density p m is lowered by larger £,’ SUI f, steep rise of e n at p >0.12 tin 3 and larger negative values of 
Pp 0) at p > 0.1 fm 3 compared to various Skyrme-Hartree-Fock calculations. 

3. For supernova matter (x ~ 0.3) 

Relatively low e n at low neutron densities <0.1 fm -3 acts to favor mixed phases rather than the uniform phase 
and Zisurf acts in the opposite way in comparison with the Skyrme-Hartree-Fock result by RBP 0. 

III. SIMULATIONS AND RESULTS 

Using the framework of QMD, we have solved the two major questions posed in the beginning of this article 
SEES . In the present section, we will review these works 1 . Hereafter, we set the Boltzmann constant = 1. 

In our simulations, we treated the system which consists of neutrons, protons, and electrons in a cubic box with 
periodic boundary condition. The system is not magnetically polarized, i.e., it contains equal numbers of protons (and 
neutrons) with spin up and spin down. The relativistic degenerate electrons which ensure charge neutrality are regarded 
as a uniform background IHl (see Refs. for effects of the electron screening). The Coulomb interaction is 

calculated by the Ewald method taking account of the Gaussian charge distribution of the proton wave packets. 


1 This section is based on our recent review article Hi. 









FIGURE 3. (Color) Nucleon distributions of the pasta phases in cold matter at x = 0.5; (a) sphere phase, O.lpo (Lbox — 
43.65 fm, N = 1372); (b) cylinder phase, 0.225p 0 (L box = 38.07 fm, N = 2048); (c) slab phase, 0.4p 0 (L box = 31.42 fm, N = 2048); 
(d) cylindrical hole phase, 0.5po (L box = 29.17 fm, N = 2048) and (e) spherical hole phase, 0.6po (L b ox = 27.45 fm, N = 2048), 
where L box is the box size and N is the total number of nucleons. The whole simulation box is shown in this figure. The red particles 
represent protons and the green ones neutrons. Taken from Ref. [39|. 


1. Realization of the Pasta Phases and Equilibrium Phase Diagrams 

In Refs. 138, 39], we have reproduced the pasta phases from hot uniform nuclear matter and discussed phase 
diagrams at zero and finite temperatures. In these works, we first prepared a uniform hot nucleon gas at the temperature 
T ~ 20 MeV as an initial condition, which is equilibrated for ~ 500 — 2000 fm !c in advance. To realize the ground 
state of matter, we then cooled it down slowly until the temperature got ~ 0.1 MeV or less for <9(10 3 — 10 4 ) fm/c, 
keeping the nucleon number density constant. In the cooling process, we mainly used the frictional relaxation method 
(equivalent to the steepest descent method), which is given by the QMD equations of motion plus small friction terms. 
In the case of finite temperatures, we also used thermostat to reproduce the equilibrium states. 

The resultant typical nucleon distributions of cold matter at subnuclear densities are shown in Fig. [3] for proton 
fraction of matter x = 0.5. We can see from these figures that the phases with rodlike and slablike nuclei, cylindrical 
and spherical bubbles, in addition to the phase with spherical nuclei are reproduced. The above simulations have shown 
that the pasta phases can be formed dynamically from hot uniform matter within a time scale of ~ (9(10 3 — 10 4 ) fm/c. 

We show snapshots of nucleon distributions at T = 1,2 and 3 MeV for a density p = 0.225po in Fig.0 This density 
corresponds to the phase with rodlike nuclei at T = 0. From these figures, we can see the following qualitative features; 
at T ~ 1.5—2 MeV (but snapshots for T ~ 1.5 MeV are not shown), the number of the evaporated nucleons becomes 
significant; at T > 3 MeV, nuclei almost melt and the spatial distribution of the nucleons are rather smoothed out. 

When we try to classify the nuclear structure systematically, the integral mean curvature and the Euler characteristic 
(see, e.g.. Ref. fiTl and references therein) are useful. Suppose there is a set of regions R, where the density is higher 
than a threshold density p tb . The integral mean curvature and the Euler characteristic for the surface of this region 
dR are defined as surface integrals of the mean curvature H = + Kb)/2 and the Gaussian curvature G = K\ Kb, 

respectively; i.e., J dR HdA and % = ^ [ dR GdA, where K'| and Kn are the principal curvatures and dA is the area 
element of the surface of R. The Euler characteristic % depends only on the topology of R and is expressed as % = 
(number of isolated regions) — (number of tunnels) + (number of cavities). Using a combination of these two quantities 




(a) T=1 MeV 



(b) T=2MeV 



(c) T=3MeV 



FIGURE 4. (Color) Nucleon distributions for* = 0.5, p = 0.225po at the temperatures of 1, 2 and 3MeV. The total number of 
nucleons N = 16384 and the box size L^ox = 76.14 fm. The upper panels show the top views along the axis of the cylindrical nuclei 
at T = 0, the lower ones the side views. Protons are represented by the red particles, and neutrons by the green ones. Taken from 
Ref. m. 



FIGURE 5. (Color) Phase diagram of matter at x = 0.5 plotted in the p - T plane. The dashed and the dotted lines on the diagram 
show the phase separation line and the limit below which the nuclear surface can be identified, respectively. The dash-dotted lines are 
the phase boundaries between the different nuclear shapes. The symbols SP, C, S, CH, SH, U stand for nuclear shapes, i.e., sphere, 
cylinder, slab, cylindrical hole, spherical hole and uniform, respectively. The parentheses (A,B) show intermediate phases between 
A and B-phases with negative %, which are different from coexistence phases of them. The regions (a)-(g) correspond to the nuclear 
shapes characterized by Jg R HdA and X as follows: (a) J^ R HdA > 0, * > 0; (b) J dR HdA > 0, X = 0; (c) J dR HdA > 0, X < 0; (d) 
J dR HdA = 0, x = 0; (e) J dR HdA < 0, X < 0; (U L R HdA <0, X =Q- (g) J dR HdA < 0, X > 0. Simulations have been carried out 
at points denoted by circles. Adapted from Ref. |38] 














































calculated for nuclear surface 2 , each pasta phase can be represented uniquely, i.e., for the phase with spherical nuclei: 
JdR HdA > 0, X > 0, cylindrical nuclei: J dR HdA >0, X = 0, slablike nuclei: J dR HdA = 0, X = 0. cylindrical bubbles: 
fdR HdA < 0, X = 0, and spherical bubbles: f dR HdA < 0, X > 0. We note that the value of % for the ideal pasta phases 
is zero except for the phase with spherical nuclei or spherical bubbles with positive negative % is not obtained for 
the pasta phases. 

The phase diagram obtained for x = 0.5 is plotted in Fig. 0 As shown above, nuclear surface can be identified 
typically at T < 3 MeV (see the dotted lines) in the density range of interest. Thus the regions between the dotted line 
and the dashed line correspond to some non-uniform phase, which is however difficult to be classified into specific 
phases because the nuclear surface cannot be identified well. 

In the region below the dotted lines, where we can identify the nuclear surface, we have obtained the pasta phases 
with spherical nuclei [region (a)], rodlike nuclei [region (b)], slablike nuclei [region (d)], cylindrical holes [region 
(f)] and spherical holes [region (g)]. It is noted that in addition to these pasta phases, structures with negative % have 
been also obtained in the regions of (c) and (e); matter consists of multiply connected nuclear and bubble regions (i.e., 
spongelike structure) with branching rodlike nuclei, perforated slabs and branching bubbles, etc. A detailed discussion 
on the phase diagrams is given in Ref. & 


2. Structural Transitions between the Pasta Phases 

In Ref. |H , we have approached the second question asked at the beginning of this article. We have performed QMD 
simulations of the compression of dense matter and have succeeded in simulating the transitions between rodlike and 
slablike nuclei and between slablike nuclei and cylindrical bubbles. 

The initial conditions of the simulations are samples of the columnar phase (p = 0.225po) and of the laminar phase 
(p = 0.4po) of 16384-nucleon system at x = 0.5 and T — 1 MeV. These are obtained in Ref. |39j], which are presented 
in the last section. We then adiabatically compressed the above samples by increasing the density at the average rate 
of ~1.3 xlO 5 po/(fm/c) for the initial condition of the columnar phase and —7.1x10 6 po/(fm/c) for that of the 
laminar one. According to the typical value of the density difference between each pasta phase, — O.lpo (see Fig. 
0 , we increased the density to the value corresponding to the next pasta phase taking the order of 10 4 fm/c, which 
was much longer than the typical time scale of the nuclear fission, — 1000 fm/c. Thus the above rates ensured the 
adiabaticity of the simulated compression process with respect to the change of nuclear structure. Finally, we relaxed 
the compressed sample at p = 0.405po for the former case and at 0.490po for the latter one. These final densities are 
those of the phase with slablike nuclei and cylindrical bubbles, respectively, in the equilibrium phase diagram at T — 1 
MeV (see Fig. |5j). 

The resulting time evolution of the nucleon distribution is shown in Figs.|6]andQ As can be seen from Fig. |6] the 
phase with slablike nuclei is finally formed [Fig.|6j-(6)] from the phase with rodlike nuclei [Fig.|6j-(1)]. The temperature 
in the final state is — 1.35 MeV. It is noted that the transition is triggered by thermal fluctuation, not by the fission 
instability: when the internuclear spacing becomes small enough and once some pair of neighboring rodlike nuclei 
touch due to thermal fluctuations, they fuse [Figs. ER2) and[U(3)]. Then such connected pairs of rodlike nuclei further 
touch and fuse with neighboring nuclei in the same lattice plane like a chain reaction [Fig.|6j-(4)]; the time scale of the 
each fusion process is of order 10 2 fm/c, which is much smaller than that of the density change. 

The transition from the phase with slablike nuclei to the phase with cylindrical holes is shown in Fig. Q When 
the internuclear spacing decreases enough, neighboring slablike nuclei touch due to the thermal fluctuation as in the 
above case. Once nuclei begin to touch [Fig. □(2)] bridges between the slabs are formed at many places on a time 
scale (of order 10 2 fm/c) much shorter than that of the compression. After that the bridges cross the slabs nearly 
orthogonally for a while [Fig.0-(3)]. Nucleons in the slabs continuously flow into the bridges, which become wider 
and merge together to form cylindrical holes. Afterwards, the connecting regions consisting of the merged bridges 
move gradually, and the cylindrical holes relax to form a triangular lattice [Fig. 13(6)]. The final temperature in this 
case is — 1.3 MeV. 

Trajectories of the above processes on the plane of the integral mean curvature J dR HdA and the Euler characteristic 
X are plotted in Fig. [8] This figure shows that the above transitions proceed through a transient state with “spongelike” 
structure, which gives negative X- As can be seen from Fig. [3(a) [Fig. |8}(b)], the value of the Euler characteristic 


Nuclear surface generally corresponds to an isodensity surface for the threshold density p t h ~ 0.5po in our simulations. 




FIGURE 6. (Color) Snapshots of the transition process from the phase with rodlike nuclei to the phase with slablike nuclei (the 
whole simulation box is shown). The red particles show protons and the green ones neutrons. After neighboring nuclei touch as 
shown by the circle in Fig.[6j(2), the “compound nucleus” elongates along the arrow in Fig. Ek3). The box size is rescaled to be 
equal in this figure. Adapted from Ref. m. 


begins to decrease from zero when the rodlike [slablike] nuclei touch. It continues to decrease until all of the rodlike 
[slablike] nuclei are connected to others by small bridges at f ~ 9840 fm/c [~ 12000 fm/c]. Then the bridges merge 
to form slablike nuclei [cylindrical holes] and the value of the Euler characteristic increases towards zero. Finally, the 
system relaxes into a layered lattice of the slablike nuclei [a triangular lattice of the cylindrical holes]. Thus the whole 
transition process can be divided into the “connecting stage” and the “relaxation stage” before and after the moment 
at which the Euler characteristic is minimum; the former starts when the nuclei begin to touch and it takes ~ 3000 - 
4000 fm/c and the latter takes more than 8000 fm/c. 


3. Formation of the Pasta Phases 

In closing the present article, let us briefly show our recent results of a study on the formation process of pasta nuclei 
from spherical ones; i.e., a transition from the phase with spherical nuclei to that with rodlike nuclei. Time evolution of 
the nucleon distribution in the transition process is shown in Fig. [9] The initial condition of this simulation is a nearly 
perfect bcc unit cell with 409 nucleons (202 protons and 207 neutrons) at T ~ 1 MeV. We compressed the system in a 
similar way to that of the simulations explained in the previous section. The average rate of the density change in the 
present case is ~4.4x 10~ 6 po/(fm/c). Since the two nuclei start to touch [see the circles in Fig.|5J-(2)], the transition 
process completes within ~ 2500 fm/c and the rodlike nucleus is formed. The final state [Fig. |9]-(4)] is a triangular 
lattice of the rodlike nuclei. 

The present simulation has been performed using a rather small system; effects of the finite system size in this 
simulation should be examined. Detailed investigation of the transition using a larger system will be presented in a 
future publication [[34;]. 














FIGURE 7. (Color) The same as Fig. 1 for the transition from the phase with slablike nuclei to the phase with cylindrical holes 
(the box size is not rescaled in this figure). After the slablike nuclei begin to touch [see the circle in Fig. 0(2)], the bridges first 
crosses them almost orthogonally as shown by the arrows in Fig. 0(3). Then the cylindrical holes are formed and they relax into a 
triangular lattice, as shown by the arrows in Fig. 0(6). Adapted from Ref. [3J|. 
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FIGURE 8. (Color) Time evolution of Jg R HdA and / during the simulations. The data points and the error bars show, 
respectively, the mean values and the standard deviations in the range of the threshold density p t j, = 0.3 - 0.5po, which includes 
typical values for the nuclear surface. The panel (a) is for the transition from cylindrical (C) to slablike nuclei (S) and the panel (b) 
for the transition from slablike nuclei to cylindrical holes (CH). Transient states are shown as [C,S] and [S,CH] for each transition. 
Adapted from Ref. & 


0 

-5 
-10 
-15 
-20 

0 1000 2000 3000 4000 5000 -2000 -1500 -1000 -500 0 


(a) p 0,225-0.405p 0 4 52 gof m /c- 

/ 6050fm/c - 


A.. 


17720fm/c 

14050fm/c 


[C,S] 


i 7560fm/c 




11380fm/c 


9840fm/c 




3310fm/c 


(b) P = 0.400-0.490p 0 ofm/c> 

Q|-| 8460fm/c/ S 

19610fm/c [S,CH] 

12000fm/c 

_..Jh 9160fm/c 


27370fm/c \ 


’t" . * 

13510fm/c 


9870fra/c 

10590fm/c 


IV. CONCLUSION 

We approached the two questions posed in Section I using the framework of QMD. According to the results of our 
simulations, our answer is strongly affirmative for both questions. The nuclear interaction used in our simulations 
shows generally reasonable properties at subnuclear densities not only for symmetric nuclear matter but also for 
neutron matter. This result also supports our conclusion. 
















FIGURE 9. (Color) Snapshots of the transition process from the bcc lattice of spherical nuclei to the triangular lattice of rodlike 
nuclei (the whole simulation box is shown). The red particles show protons and the green ones neutrons. The box size is rescaled 
to be equal in this figure. 
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